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ABSTRACT 

o 
o 

CN ' We consider element diffusion during the hnear growth of structure in the intergalac- 

tic medium (IGM) before the redshift of reionization. The elements produced during 
• big bang nucleosynthesis (such as D, ^He or "^Li) condensed in collapsing regions more 

(<~^ , readily than hydrogen. Diffusion was most effective when the IGM was neutral since 

, the coUisional cross-section between neutral species is smaller than the corresponding 

Coulomb cross-section in an ionized gas. We find that diffusion led to small deviations 
I in element abundances of up to 0.1%. Abundance measurements to this level of preci- 

' sion in low-metallicity environments could provide information about the thermal and 

' ionization history of the primordial gas during the early formation of structure in the 

^ ' IGM. 



i-^l Subject headings: cosmological parameters — cosmology: theory — diffusion — galax- 

I ' ies: abundances 

o 



1. Introduction 

^ ' Recent attempts to measure the primordial ^He abundance report accuracies better than ~ 1% 

■ (Izotov & Thuan 1998; Olive, Steigman, & Skillman 1997; Peimbert, Peimbert, & Ruiz 2000). Such 

accurate measurements are useful for testing the consistency of big bang nucleosynthesis and also 
for estimating the mean baryon density, 17;, (Schramm & Turner 1998). The primordial ^He mass 
fraction, Y, is often inferred from the emission lines of blue compact galaxies (e.g., Peimbert & 
Torres-Peimbert 1974; Pagel, Terlevich, &; Melnick 1986) under the implicit assumption that it 
had the same initial value throughout the Universe, including within underdense regions of the 
intergalactic medium (IGM). 

However, a net difference between the element abundances in overdense and underdense regions 
could arise if diffusion occurs during the gravitational collapse of massive objects. Different atomic 
species, because of their different masses, will have different thermal speeds and will therefore 
experience different pressure forces per unit mass. Consequently, different elements would not 
flow into an overdense region at the same rate; instead, heavier elements would obtain higher 
relative abundances in galaxies than in voids. Diffusion is generally moderated by a frictional force 
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arising from inter-species collisions. These collisions are most frequent for an ionized gas (Spitzer 
1956) in which the characteristic cross-section for ion-proton collisions, ~ 'ir{Ze'^ /kBT)'^ In Aq ~ 
2 X 10^^^Z(T/10^ K)~^ cm^, is far greater than the collision cross-section for the corresponding 
atoms in a neutral gas, ~ 10^^^ cm^ (here, Ze is the ion charge, T is the temperature, ks is 
Boltzmann's constant, and In Ac ~ 20 is the Coulomb Logarithm). We therefore focus our 
attention on the period between cosmological recombination and reionization during which the 
IGM was neutral (Loeb & Barkana 2001). During this stage, atomic collision cross sections may 
be small enough to allow some diffusion to occur. 

In this paper we quantify the diffusion of ^He, ^Li, and D relative to that occurs as a spher- 
ical, slightly overdense region begins to collapse. We assume a ACDM cosmology, characterized 
by densities relative to critical of 0,^ = 0.3, f^A = 1 — ^l^, and Jlfe = 0.045 for matter, vacuum, 
and baryons; a Hubble constant of Hq = 70 km Mpc~^; a ^He primordial mass fraction of 
Y = 0.244; and primordial ^Li and D number densities relative to of Xl = 1-75 x 10^^*^ and 
Xf) = 2.0 X 10^^, respectively. The effects of diffusion are most important where the gas pressure 
matters, namely around the comoving Jeans length scale Aj = {2ii/kj). For redshifts z ^ 140, 
when the baryonic temperature is primarily determined by the adiabatic expansion, the Jeans 
wavenumber is (Barkana & Loeb 2001) 

21 1/2 



kj 



0.9 (1 + z) fXTTlH^rnHf 

ksT (z) 



(1) 



where h = {Hq/WO km Mpc~^), fj, is the mean molecular weight, ttih is the mass of hydrogen 
and T (z) the baryonic temperature at a redshift z. 

In § ^, we derive the differential equations used for our analysis. These are linearized equations 
which are strictly valid only as long as the overdensities remain small. However, we assume that 
reionization occurs before the collapsing region becomes significantly nonlinear and we stop our 
calculation at reionization because charged particle cross sections suppress any further diffusion. 
In § 1^ we discuss the initial conditions used for our differential equations. In § ^, we compute both 
the mean abundance deviation and the typical abundance fluctuation within a collapsing region. 
We then use a Press-Schechter mass function to average the resulting deviations over all collapsed 
objects in the Universe. In § |5|, we compare our model to recent measurements of ^He abundances. 
Finally, we present our conclusions in § 



2. Basic Equations 

We first establish our notation and illustrate our assumptions by carrying out a brief derivation 
of the equations satisfied by density perturbations. Each particle species has a mass density, ps, 
and a bulk velocity, Vg , which are related by the continuity equation, 

dps 
dt 



+ Vp • (psVs) = 0, (2) 
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where we start by using physical (proper) coordinates. Since the Universe is initiahy homogeneous, 
Ps = Ps (t) and all particles move with the Hubble flow: = H{t)x (Peebles 1993). In this 
approximation, equation (0) becomes 



— + ^p.H = 0. 



(3) 



We next allow for density and velocity perturbations through the definitions 

Ps{^,t)=ps{t)[l + 6s{^,t)], (4) 

d'vi (X, t) = a^J (X, t) + ioi^ (X, t) + [H it) + 5H, (x, t)] Cnecfcer- (5) 

The velocity gradient is decomposed into a symmetric, traceless tensor (u*-'), an antisymmetric ten- 
sor {u>l^), and a trace. We insert these definitions into equation (§) and linearize the perturbations 
to obtain 

-^ + {^r,.V)5s + 3SHs=0. (6) 
ot 

Switching from now on from physical to comoving coordinates, and letting an overdot represent a 
time derivative, equation (^) may be compactly expressed as 



-36 Hs. 



(7) 



The density and bulk velocity also satisfy a force equation; in comoving coordinates, 

v, = p;^ ^ F,. 

forces 



(8) 



The sum on the right-hand side includes all forces acting on species s. We linearize this equation 
and make use of equation (0) to obtain a new equation which does not explicitly contain the velocity 
perturbations. 



5s + 2Hds - 3— - 3^2 
ot 
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p 



forces 



(9) 



We now switch independent variables from time to redshift via the Friedmann equation, which 
leads to 

2 



Hot 



: sinh 



-1 



(1 + ^) 



-3/2 



(10) 



However Q\ = 0.7 only becomes significant for z ^ 1; in this paper we consider redshifts z ^ I for 
which 



Hat 



3^A^r 



:i + z) 
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Denoting a derivative with respect to z by a prime, equation (P) becomes 
1 3 / dH 1 



S' + 



2(l + z) 



Qn^H^ [l + z 



flmH^ (1 + zf 



V • 



' forces 



(11) 



(12) 
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We now obtain expressions for the forces. The force due to gravity, when expressed in terms 
of the potential, (/>, is = —psVcj) . We insert this into equation (12) and use Poisson's equation 
which gives 



2{l + z) ' a^H^{l + zr 



-E', fJ + ^ (13) 

Ps ^forces 2Vlrn ^ 

species 



where the primed sum means a sum over all forces except gravity. 

Although all types of particles are affected by gravity, we assume that only the baryons ever 
obtain a non-zero temperature. Each species of baryons experiences a different pressure force which 
is determined by the negative of its pressure gradient, 

F, = - VP, = - — {psVTs + T.Vps) , (14) 
rus 

where we have used for the atomic mass and Tg for the temperature of species s. Dividing 
equation ( p!4[ ) by the mass density and taking the divergence, we see that 

y.l± = .h^{y^Ts + T,VH,) (15) 
Ps rus^ / 

to first order in the perturbations. We define Tg to be the temperature in the absence of perturba- 
tions, and express the effects of perturbations in terms of the small quantity r,: Tg = Ts{l + Tg). 
Inserting this into equation (|l5|) and linearizing once more we see that 



V- — 



Fs kefg 



Ps rng 



(VV, + V^6g) . (16) 



The temperature of the baryons is determined by two effects, the coupling of free electrons to 
the background radiation and the adiabatic expansion (Peebles 1993). This gives 

^T, ^ ^Sj^ 2T.a^^ (17) 

at 1 + X SnieC cips at 

where x is the free electron fraction, or is the Thomson cross section, a is the Stefan-Boltzmann 
constant, me is the electron mass, c is the speed of light, and is the temperature of the cosmic 
microwave background (CMB). Switching the independent variable to z and linearizing, equation 
( [TTD becomes 

dTg^ ^ X SaraT^jl + zf^ / _ Tg \ 2Tg 2^ ^, 

dz l + x 2,mecHQ./rL;^ \ Tq{1 + z)) (1 + z) 3 ^ ' 

where Tq = 2.73 K is the CMB temperature at present. 

The mean baryon temperature, T,, is defined to be a solution of 

X SaraT^ (1 + zf^ / ^ I (19) 

dz l + x 3meCHQ^yQrn \ Tq {1 + z) I {1 + z) 



- 5 - 



This indicates that Ts is independent of baryonic species and so all baryons would be at the same 
temperature in the absence of perturbations. We use the computer program RECFAST (Seager, 
Sasselov, & Scott 1999) to determine the free-electron fraction as a function of redshift. Then we 
use equation ( p!9[ ) to compute the mean baryon temperature. Both quantities are shown in Figure |l|. 

After linearizing equation (|l8|), we obtain the equations for the temperature perturbations, 

X SaraT^ (1 + zf^ Ws ^ x IQaraTl (1 + zf'^ 6s ^^q) 
1 + X SiriecHo^/Qm Ts 1 + X 9'meCHo\/flm Ts^ 

where Ws = Ts — ^6s- With this definition of Wg, equation (|l^) becomes 

V.!^ = -MfV<5. + wV (21) 



Ps rus 

Because of the spatial derivatives, we are forced to work in comoving Fourier space. However, for 
simplicity of notation we will continue to simply write 6s and Ws, although from now on we use 
these symbols to represent the Fourier components of the perturbations. 

Finally we take into account the frictional term arising from inter-species collisions. Physically, 
this force is non-zero only if two species have different bulk velocities; thus, the force on one species 
due to a second species may be written as Fi = k (x, t) (v2 — vi) (Burgers 1969). Since the 
difference in bulk velocities is already linear in the perturbations, we only need to calculate the 
coefficient k to zeroth order for the linearized equations of motion. This allows us to calculate 
K under the assumption that all baryons share the same temperature. In general, the frictional 
force per unit volume acting on species i due to species j is obtained by multiplying the number 
of collisions per unit time per unit volume by the momentum transfer per collision. We multiply 
this by the probability that species i and j have initial momentum pi and pj and that the center 
of mass scattering angle is x- Finally, we integrate over all parameters, 

1 dai 



Fij = J d^p, J d-^pj J dn niUjaij \vi - Vj\Api (pi, Pj, x) fi (Pi) fj (Pj) (?»' Pi' • (^2) 



We assume that the collisions are elastic and that the atoms can be approximated by hard 
spheres. We take each species to have a Maxwellian velocity distribution. Following the integration 
equation ( p2|) becomes 



IGrriij ksT f°° = _~2 / kBT\ 
Fij = - -— i W ninjaci{-Vi-Vj) dxx e Qd\x\ . (23) 



3V27r ^ rriij ' " Jo 

Here, rriij = rriimj / {rrii + rrij) is the reduced mass, is the classical hard-sphere cross section, 
iT-i = pi/rrii, and the integral is a quantum mechanical correction (Massey, Burhop, & Gilbody 
1969). The function (q) is given by 

^ oo 

= — ^ [4 {21 + 1) sin2 r]i - 2 {I + cos 2r]i - cos 2r]i+i + cos 2 {7]i - r?;+i))] . (24) 
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The integration variable x is related to q hy q = xay 2mijkBT /fi^ , and the phase shifts are given 

by Vo = —Q and tanry^ = (— '/i+1/2 ('i')/'^-«-i/2 (^) > where the functions Jn{q) are Bessel 
functions of order n. The distance scale a represents the sum of the radii of two colliding atoms; 
we adopt the values rn = = 0.53A , rne = 0.3lA , and rn = 1.45A (Clementi, Raimondi, 
&: Reinhardt 1963; Slater 1964). We evaluate the integral in equation (|2^ ) numerically and plot 
the results in Figure § for each type of inter-species collision. Quantum mechanics becomes more 
important as the temperature decreases but the correction is never greater than a factor of four. 
For brevity, we henceforth denote this integral by lij (z). 



We now summarize the equations governing each species. For dark matter equation (12) 

s'L + ^ (1 + = I (1 + E (25) 

species 

For each baryonic species there are also pressure and collisional friction terms, and so 



5^ + ^{i + zr'6i = Ui + z)-' E ?r^. (26) 



j=all species 

(1 + z)"^ k'^ ( ^6i + w, 



kpT 3 
m.f]^//2^ V3' 



'2^V ^^"'^ 



\Jmimj {mi + mj) 



where /c is the comoving wave number. From equation (26) we can easily find the differential 
equations satisfied by the differences relative to hydrogen, = 5i — 5h'- 



Ar + ^(l + .)-^A^ = 7^(1 + ^)-^^^ 



5 f 6h ^ Aj -6h \ ^{wh Wi 



3 \mH rrii J ymn fn. 



(27) 



""G V irT3 ^ / , ^ ^ (1 + (^) ^ 

V iim" ^yj:^ JuiHinj [mH + rrij) 



3. Initial Conditions 

We consider dark matter perturbations initially of the spherical top-hat form (x, ^im) = 
6ini9 {R — |x|), which has the Fourier components 

471" (5', • 

5dm (k, Zini) = ^3"' (sin kR - kR cos kR) . (28) 

Here, and R are constants and represents the step function. We take the initial values of 
A He, An, Ad and their derivatives to be zero, but we would like to choose the initial values of 
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^dm' ^Hi and that correspond to the most rapidly growing mode. To gain some insight into how 
to make this choice, we consider a simpHfied problem where baryonic diffusion is negligible. We 
approximate the baryonic temperature by the temperature of the background radiation (valid for 
z ^ 200; see Figure |l]). Equation (26) for the baryon overdensity becomes 



dm 



(29) 



We use this equation to express 6drm differentiate, and insert the results into equation (p5|). The 
outcome is a fourth order differential equation for 6h. It has four independent, power law solutions 
which we write as 6^ = J2i=i (1 + -z)^'- The four Ai are the four arbitrary constants associated 
with a linear fourth order equation. The four pi are given by 



Pi 




2 



fc£ I 1 
2 2 



(30) 



A;2 



2 



where is the squared comoving wavenumber normalized by {fiQmHQ /kBTo). The resulting 
expression for the dark matter overdensity is 



with the constants Cj given by 



Pi 



(31) 



i=l 



( 2Vir 



\ 3r2c(^ 



J. 



\ 3f^dri 



+ 



(32) 



Motivated by this solution, we choose our initial conditions so that the perturbation will 



grow like the most negative (i.e., the third one listed in equation 30). We let the three Ai 
corresponding to all but the largest growth rate be zero, and let the final Ai be determined by the 
tophat distribution (equation [2^) . The remaining initial conditions may now be written in terms 
of the known Cj and pj. 



4. Calculations 
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4.1. Evolution of the k-Modes 

We define a new set of variables, 

^3 

^ Andini {sin kR — kR cos kR) ^' ^ ^ 

" AirSini {sin kR-kR cos kR) ~ ' ^^^^ 

2^ — rJJJ ( OO 1 

AirSini {sin kR — kR cos kR) ^ 
Tliese quantities are special because they not only satisfy the same differential equations as the 
5 and w, but also satisfy simpler initial conditions, 6dm = 1 > ^'dm ~ i+z- > = ^ , and 
K = — — V • The parameters (Sini,R) have disappeared from the problem. Therefore we only 
evolve each normalized /c-mode once; we obtain the physical overdensity for a given 6ini and R by 
multiplying each /c-mode by the appropriate factor from equation (|3^). 

We have set Zini = 10^ following cosmological recombination. The final redshift is marked by 
the first of two possible events. Our equations are applicable in the linear regime, so we must stop if 
the object becomes nonlinear. At a transition to nonlinearity, the density would be greatly enhanced 
which would result in the suppression of diffusion. Therefore, we expect diffusion to be effective 
in the linear regime where our equations are valid. The second stopping condition relates to the 
reionization of the Universe. At that time, large Coulomb cross sections between charged particles 
are able to effectively halt diffusion. In conclusion, we consider the case where the object becomes 
nonlinear after reionization and adopt Zreion — 

8 (Haiman & Loeb 1998, 1999a,b; Gnedin & Ostriker 
1997; Chiu & Ostriker 2000; Gnedin 2000) as the stopping point of our diffusion calculation. Since 
the collapsed fraction of the Universe is small at z = 8 (Barkana & Loeb 2001), most of the material 
within galaxies that form much later (at ^ 3) was still in the linear regime at z 8, where our 
diffusion calculation is adequate. 

Several results of having carried out this prescription are shown in Figure ^. All length scales 
have been measured in terms of kj^ , defined in equation (|l]). This expression involves the baryonic 
temperature which is set by the adiabatic expansion for redshifts less than zt ~ 138, T ( 
To {I- 
value 



To (1 + z)^ (1 + Zt) ^; as discussed in § 0. At Zreion, the Jeans wavenumber therefore obtains the 



kj = 1040 Mpc" 



9 \ / nmh^\ /2.73 K 



l + zj \ 0.147 / V Tf 







1/2 

(36) 



which is the value used in Figure ^. We see that the dark matter overdensity has evolved approx- 
imately proportionally to (1 -|- z)~^ , as might have been expected since gravity is the only force 
applicable to dark matter. On large scales corresponding to k ^ kj, the baryons follow the dark 
matter; however, on small scales corresponding to k ^ kj, pressure forces prevent baryonic modes 
from collapsing. This leads to the slight decline in dark matter amplification near the Jeans scale 
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and the dramatic decline in the ampUfication of baryonic modes there. In the lower panel of Figure 
3, we see that the perturbations of the different baryonic species do not experience exactly the 
same amplification. At <C fcj, pressure is negligible for all species, and the modes of each are 
amplified. Since hydrogen is the lightest species, it reaches its Jeans scale at the lowest k. Thus, 
all of the A^^ become positive at approximately the same wavenumber. As k increases, we begin to 
reach the Jeans scales of the heavier species. Then the growth of A;-modes of this second species 
is also suppressed, and the A^; approach zero as k gets large. However, Figure ^ indicates that 
the overall magnitude of diffusion is not solely determined by the Jeans scale of a species, or else 
the effect would have been largest for ^Li. This is because of the frictional term which serves to 
suppress diffusion, which is proportional to the square of the radius of an atom and therefore the 
greatest suppression is associated with the largest atoms. 

The normalized overdensities in Figure ^ were obtained by evolving 1025 A;-modes, equally 
spaced between with < k/kj < 100. To obtain the physical overdensities, each mode must be 
multiplied by its normalization factor. The normalization factor depends on the object assembly 
radius R, and is fairly rapidly varying compared to the normalized overdensities. To get an accurate 
determination of the physical overdensity, we interpolate the normalized overdensity so that we 
sample it at a rate larger than the frequency of the normalization factor. We then use equation 
(|3^ to obtain the physical fe-modes. One such sampling is shown in Figure The value of R 
used for this plot corresponds to an enclosed mass of M = 10^ M0 = pcrit^m^TTR^. The relative 
errors incurred in interpolation are typically ~ 10~^, and seldom exceed 10~^ with the 1025 initial 
points. The normalization factor is also proportional to the initial dark matter perturbation, 6ini, 
which we express in terms of the collapse redshift of the object (zcou) through the approximate 
relation 6ini (1 + Zini) = 5coU (1 + Zcoii) ■ We take 5coU = 1-69 for the overdensity of collapse and, 
unless otherwise stated, ZcoU = 2.5 for the redshift of collapse. Such objects would be just on the 
fringe of nonlinearity at Zreion = 8. Since diffusion within all objects was shut off at Zreion, the 
effects of diffusion would be decreased for objects that collapsed later. 



4.2. Predicted Statistics 

4-2.1. Mean Abundances 

Next we quantify the extent to which element abundances in collapsed objects differ from their 
primordial values. Abundances of ^He are commonly reported in terms of the mass fraction, 

PHe (1 + ^He) + pH {'^ + 5h) ' 

where we have neglected the contributions of all species but ^He and to the mass of the Universe. 
After linearizing this equation, we see that 

Y -Y 



Y 



(1 - Y) {6He - 6h) ■ (38) 
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in contrast, the abundances of ^Li and D are typically reported in terms of their number densities 
relative to hydrogen. In this case, 

"""'^ " ^"'^ = {6u,D - 5h) . (39) 



Xl. 



D 



To see how the element abundances in collapsed objects differ, on average, from their primordial 
values, we average equations (p8|)- (|3^) over the original sphere of comoving radius R, 



= -1- / d^x / ^Ake*'^''. (40) 

iTlR^ J sphere J (lirf ^ ' 



After integrating, we see that 



3 dk 

Ar = — Afe (sin kR - kR cos kR) . (41) 



We carry out this final integration numerically. As described in § 4.1, we evolve 1025 /c-modes 
between and WOkj to get the A^; (where x could denote any of ^He, ^Li, D). We interpolate to 
get values for A^^ at a fixed number of points, and then multiply each point by its normalization. 
We then evaluate the integral in equation (^Tj) with these points. Then we repeat the procedure, 
but with twice as many interpolated points. If the two results differ by less than a fraction of 
5 X 10^^, we stop and accept the result. 

Figure |^ illustrates the dependence of these deviations on the size of the collapsed object. 
Instead of measuring the size of the object by its assembly radius, R, we instead measure it in 
terms of its mass defined by 

M = Pcrit^m^7rR\ (42) 

We see that the deviation from the primordial mean shrinks as the collapsed object increases in 
size approximately according to a power law with slope —0.33. A simple argument yields the same 
result: assuming that the collapsing region starts out with primordial abundances throughout its 
volume and that diffusion occurs through the boundary of the region, the deviations from primordial 
abundances should fall off like 1/R, or M~^/^. 

Now that element abundances for particular collapsed objects are known, we proceed by aver- 
aging these abundances over all objects in the Universe that have collapsed by a particular redshift. 
This statistic is approximated by using a Press-Schechter mass function; for ''He, 

/?? YmM^cIM 

M^dM ' ^ ^ 



where dn/dM is the number density of objects of mass M that have collapsed at a given redshift 
and Ym is the averaged ^He mass fraction within an object of mass M. Equation (42) is also 



appropriate for the ^Li and D abundances if we replace Y with X^^d. Note also that there is 



a 
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minimum mass scale of galaxies in which star formation is possible that we denote by Mmin (Abel 
1995; Tegmark et al. 1997). We use a code developed by Eisenstein & Hu (1999) and Barkana 
(2000) to calculate dn/dM at a redshift of 2.5. The averaged abundances differ from primordial 
by the amounts shown in Figure ^. We also consider a collapse redshift of ZcoU = 0; in this case 
Figure ^ shows that the abundance deviations decrease somewhat. This is not only due to the 
resulting decrease in 5ini, but also due to the dependence of dn/dM on collapse redshift. At higher 
z, dn/dM tends to favor lower mass scales where the abundance deviations are larger (Figure ||). 



4-2.2. Fluctuations of Abundances Within Objects 

The element abundances of collapsed objects can differ from primordial abundances in two 
ways. First, the abundances of an object could be uniformly different from the primordial values. 
This possibility was already discussed in the previous section. Second, if the element abundances 
within an object vary greatly from point to point, differences from the primordial abundances may 
be observed because of the fluctuations. We consider this latter effect by computing the standard 
deviation of the abundances within a given object. The lowest-order formula for this quantity is 

a'y = Y^l-Yfai^ (44) 

for ^He and 

4 = (45) 
for 'Li and D. We have already calculated the mean values of the abundances, so now we only need 
the mean squares. We proceed as above, but this time we obtain a double integral that must be 
evaluated numerically, 

(A^) ikkA, r mA, f"" '* -," ^ - ") . (46) 

V Jm Stt'^R^ Jo Jo \ k-l k + l J ^ ^ 

The integration region is the first quadrant of the k — I plane. Since the integrand is symmetric 
with respect to the interchange of k and I, it is equal to the integral 

3 f°°rr,^ /■'' n,A f siu (k - I) R siii (k + I) R\ 

' dkkAk / dllAi ) ; , / . (47) 



-2 



M iTT^RS Jo Jo \ k-l k + l 



We carry out the integrations as with equation ( pl| ) and our results are plotted as a function 
of mass scale in Figure ^. Comparing Figures |5| and |^, we see that the relative importance of 
abundance fluctuations depends on the mass of the object in question. On small mass scales, 
abundances differ from their primordial values by an approximately uniform amount; but as the 
object size increases fluctuations from the mean are more likely to produce a difference from the 
primordial values. We also average the standard deviation over all collapsed objects in the Universe. 
Our prescription for this is identical to that given by equation (^3|), with cja/ taking the place of 
Ym ■ These results are plotted versus the minimum mass scale for galaxies in Figure Comparing 
this to Figure |6|, we again see that the relative importance of fluctuations increases with increasing 
mass scale. 
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5. 



Observational Data 



Though much recent work has focused on uncovering the primordial Li and D abundances 
(e.g., Buries & Tytler 1998a,b; Tytler et al. 1999; Webb et al. 1997; Ryan et al. 2000), the "^He 
abundance with a reported uncertainty of about 1%, is still more accurately determined. However, 
even this accuracy is at least an order of magnitude coarser than the mean effect of diffusion for 



& Thuan (1998) to determine whether or not such fluctuations are currently detectable. 

In Izotov & Thuan (1998), the helium abundance is assumed to be linear in metallicity for low 
metallicity objects. Each data point consists of a measurement of the metallicity, Zj, a measurement 
of the helium abundance Yi, and estimations for the uncertainties in each, ej and 6i. A straight line 
is fit to the data, the intercept of which reveals the primordial (zero-metallicity) helium abundance. 
The model we use is slightly different. Because of diffusion we cannot directly associate ^He 
abundance measurements of zero-metallicity objects with the primordial '^He abundance. Not only 
do we expect an object of a particular mass to have a ^He abundance differing from the primordial 
value as shown in Figure but also we expect that different measurements of the single object 
would display an intrinsic scatter about this mean as shown in Figure |^. Therefore, we fit the data 
to the following model. If an object has metallicity Z, the probability that it will have a ^He mass 
fraction Y is denoted g {Y\Z) dY. We account for measurement uncertainties by defining h (Yi\Y) 
to be the probability of measuring Yi if the object actually has abundance Y, and / {Z\Zi) to be the 
probability that an object has metallicity Z if Zj is what is observed. Thus, the total probability 
of a measurement simultaneously obtaining Yi and Zi is 



Both Z and Y are intrinsically positive, but we here assume that the uncertainties in the measur- 
ments are small enough so that the probability distributions / and h may be taken to be Gaussians, 



low-mass objects. But, as discussed in § ^ the "^He abundance is also expected to vary from point to 
point within an object. In this section, we re-examine the ^He abundance data presented in Izotov 




(48) 




21 



(49) 
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As for our model, we assume that 



1 



{Y -mZ - h) 
2^2 
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9{Y\Z) 



V27rCT2 



exp 



(51) 



The parameters to be determined by the data are m, which contains the physics of the metallicity 
dependence; 6, which is the mean value for the ^He abundance; and cr, which represents the intrinsic 
scatter about the mean. 
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We use the integration in equation (^) to obtain 

1 1 



P{Yi\Zi 



fj2 _|_ ^2g2 _|_ r2 



: exp 



{Yi — mZi — h) 



(52) 



If we determine the metalhcity and hehum abundance for N objects, the probabihty for the entire 
measurement to occur is 

N 



N 



(53) 



To maximize this probabihty with respect to the parameters m, 6, and a, we minimize the x^.^jj^g 
function 



X 



N 

E 



{Yj - mZj - h) 
0-2 +m?e'j + 5} 



(54) 



A downhill simplex method and a direction set method yielded almost identical parameters for the 
minimum m = 4.8 x 

10-3^ h = 0.244, and a = 1.52 x 10"^. We used a Monte Carlo method to 
set confidence limits on these results. Because of random errors, the observed data points Zi and Yi 
are not a unique realization of the true values of the quantities. Instead, for a given measurement 
of Zi, there are an infinite number of Yi that could have been measured. The probabilities for 
obtaining these hypothetical measurements are assumed to be given by equation ([5^). We have 
generated 10^ of these hypothetical data sets by substituting the observed values for m, 6, and a 
into equation (|5^). For each hypothetical data set, we again minimized equation ( ^4| ) to obtain a 
new triplet of parameters m, h, and a. Figure ^ shows the resulting intercepts h plotted against the 
scatters a. The uncertainty in our best-fit sigma and intercept are large, and so the observational 
situation must improve before we could settle on the true parameters. 



6. Conclusions 

We have explored the effect of diffusion during the linear growth of density perturbations in 
the neutral IGM before the redshift of reionization and have found that element abundances within 
collapsed objects differ from their primordial values. Element diffusion leads to two related effects. 
First, the mean element abundances of a collapsed object could differ from the primordial values; 
and second, diffusion could result in abundance fiuctuations within a single object. For objects with 
a total mass below 10 '''M© the former effect dominates; but for larger objects the latter effect is 
more important. The magnitude of the effect is ~ 0.005-0.2%, depending on the mass of the object 
and the element. We have then used a Press-Schechter mass function to average these distributions 
over all collapsed objects in the Universe, and again found the net effect to be < 0.1%. This effect 
is yet too small to be detected since the best observations of element abundances currently have 
accuracies of ~ 1%. 

If the effect of diffusion on element abundances will eventually be measured with sufficient 
precision in low- metalhcity systems, the results could provide important clues about early struc- 
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ture formation in the IGM. In particular, the mean deviations within objects are sensitive to the 
minimum mass of a galaxy halo. Also, the overall magnitude of diffusion is sensitive to the redshift 
of reionization since diffusion is brought to a halt at that time. And as determinations of element 
abundances get more and more accurate, diffusion will eventually have to be accounted for in order 
to make comparisons with the predictions of big bang nucleosynthesis. 

We thank Rennan Barkana for useful comments on the manuscript. This work was supported 
in part by NASA grants NAG 5-7039, 5-7768, and by NSF grants AST-9900877, AST-0071019. 
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Fig. 1. — Top panel: The free electron fraction as a function of redshift. The Universe rapidly 
approaches neutrality after z ~ 1000. Bottom panel: The average baryon temperature as a function 
of redshift. For z > 200, the free electron fraction is still large enough so that the baryons are 
coupled to the background radiation, and so the baryonic temperature decreases like (1 + z). For 
smaller redshifts, the baryon temperature is mainly determined by the adiabatic expansion, and 
drops off like {1 + z)"^. 
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Fig. 2. — The intcr-particle cross section, in units of the classical hard-sphere cross section. Solid 
line is ^H-^He, dotted is ^H-^Li, short dash is ^H-D, long dash is ^He-^Li, dot-short dash is ^He- 
D, and dot-long dash is ^Li-D. This quantum-mechanical correction factor depends on the atomic 
radii, the reduced mass, and the temperature. 
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20 40 60 

Fig. 3. — Top panel: The normalized dark matter (solid) and hydrogen (dashed) over densities. 
For small k, both are amplified. But as k increases, pressure forces inhibit baryonic modes from 
growing. Bottom panel: The normalized differences relative to Hydrogen, A^,. The solid line has 
X =^He, the dotted line has x =^Li, and the dashed line has x =D. The relatively heavy mass and 
small atomic radius of ^He make it particularly susceptible to diffusion. 
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Fig. 4. — The difference Ah^, for an overdense region with a mass of 1O^M0. This plot is obtained 
by interpolating the smoothly varying normalized overdensities, and then multiplying each mode 
by its rapidly varying normalization factor. 
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Fig. 5. — Fractional deviations from the primordial abundances for Zreion = 8. The solid line is 
the fractional deviation for ^He, the dotted line is for ^Li, and the dashed line is for D. 
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Fig. 6. — The fractional deviations averaged over all collapsed objects in the Universe, as a function 
of the minimum allowed mass scale for galaxies. We show two cases for the collapse redshift of the 
objects. For the first case with Zcoii = 2.5, the solid line corresponds to the fractional deviation for 
^Hc, the dotted line to ^Li, and the short-dashed line to D. For the second case with ZcoU = 0, the 
long-dashed line corresponds to "^He, the dot-short dashed line to ^Li, and the dot-long dashed line 
to D. 
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Fig. 7. — Standard deviations of abundances for Zreion = 8. The solid line is the standard deviation 
for ^He, the dotted Hne is for ^Li, and the short dashed hne is for D. 
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Fig. 8. — Standard deviations averaged over all collapsed objects in the Universe as a function of 
the minimum allowed mass scale. We show two cases for the collapse redshift: (i) Zcoii = 2.5, for 
which the solid line is the standard deviation for ^He, the dotted line is for ^Li, and the short-dashed 
line is for D; (ii) Zcoii = for which the long-dashed line refers to ^He, the dot-short dashed line 
refers to ^Li, and the dot-long dashed line refers to D. 
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Fig. 9. — Results from 10^ Monte Carlo data sets based on the model described in §5. The top 
panel shows all hypothetical {bhypo,o'hypo) pairs relative to the observed values {bdatajiT^data)^ while 
the bottom panel show the 6, 730 pairs most likely to occur. The extended horizontal concentration 
of points near —0.0015 occurs because of the physical constraint that all values of a be greater than 
zero. This indicates that we have not yet reached the level of accuracy required to put a lower 
bound on a. 



